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Abstract 

Using a configuration interaction approach we study statistics of the dipole 
matrix elements (El amplitudes) between the 14 lower states with J n = 
4~ and 21st to 100th even states with J = 4 in the Ce atom (1120 lines). 
We show that the distribution of the matrix elements is close to Gaussian, 
although the width of the Gaussian distribution, i.e. the root-mean-square 
matrix element, changes with the excitation energy. The corresponding line 
strengths are distributed according to the Porter-Thomas law which describes 
statistics of transition strengths between chaotic states in compound nuclei. 
We also show how to use a statistical theory to calculate mean squared values 
of the matrix elements or transition amplitudes between chaotic many-body 
states. We draw some support for our conclusions from the analysis of the 228 
experimental line strengths in Ce [J. Opt. Soc. Am. 8, 1545 (1991)], although 
direct comparison with the calculations is impeded by incompleteness of the 
experimental data. Nevertheless, the statistics observed evidence that highly 
excited many-electron states in atoms are indeed chaotic. 
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I. INTRODUCTION 



The aim of this work is to present more evidence that excitation spectra of complex 
open-shell atoms, and probably any other atom at sufficient excitation energies, display clear 
quantum chaotic features. This phenomenon is caused by strong mixing of many-electron 
excited states by the residual two-body Coulomb interaction. It manifests in particular in a 
Gaussian statistics of the El amplitudes for these states. 

Since the time of Bohr's hydrogen atom theory atoms were considered as perfectly regular 
dynamical systems. As the classical theory of chaos evolved it became apparent that highly 
excited atomic states in the Rydberg range could become chaotic if an external field is 
applied 0, as long as the underlying classical motion is chaotic. 

On the other hand, it was also due to Bohr that the notion of compound nuclei was 
introduced in physics. The behaviour of these highly excited nuclear states is essentially 
quantum-mechanical. Nevertheless, they display a number of chaotic properties. For ex- 
ample, the statistics of their energy spectra show certain universal features, and transition 
amplitudes involving compound states obey Gaussian statistics @. To describe these prop- 
erties it was suggested by Wigner that the Hamiltonian of a compound nucleus could be 
modeled by a random matrix, and different characteristics found by averaging over ensembles 
of such matrices (see reviews [§|,U). 

The first insight into quantum chaotic properties of complex atoms was given by Rosen- 
zweig and Porter |J who analyzed experimental spectra of some neutral atoms and showed 
that in heavy open-shell atoms the spectral statistics are similar to those of compound 
nuclei. That analysis was later extended and refined in |§. Of course, the study of eigen- 
values provides valuable information about the system. On the other hand, the spectral 
statistics observed in heavy open-shell atoms are similar to those of the hydrogen atom in 
a strong magnetic field 0, or even a particle in a 2-dimensional classically ergodic billiard 
||. However, the eigenstates of these quantum systems must be completely different, and it 
is clear that the eigenvalue statistics cannot really tell us much about the origin of chaotic 
behaviour, or indeed the structure of the chaotic eigenstates. 

The first inquiry into the possibility of chaos in the eigenstates of complex atoms was 
done by B. Chirikov ||. He studied configuration compositions of eigenstates of the Ce 
atom using data from the tables [ID|, and came to the conclusion that the 'eigenfunctions 
are random superpositions of some few basic states'. Inspired by that work we conducted an 
extensive numerical study of the spectra and eigenstates of complex open-shell atoms, using 
the rare-earth atom of Ce as an example fTTHT3| . This allowed us to investigate many-body 



quantum chaos in a real system. We showed that atomic excited states are in fact similar to 
nuclear compound states and developed a statistical approach for analyzing their properties. 

Unlike eigenvalues, the eigenfunctions are not observable directly. To probe the structure 
of the chaotic eigenstates one can look at the transition probabilities or matrix elements of 
some external perturbation coupling them to each other, or to regular, simple eigenstates 
(like the ground state). The matrix elements involving chaotic eigenstates must have Gaus- 
sian statistics. We showed that its main characteristics - the mean squared value of the 
matrix element between the chaotic multiparticle states (compound states), can be calcu- 
lated in terms of statistical parameters of the eigenstates and single-particle amplitudes and 
occupation numbers of the orbitals present in the compound states 
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In this work we have chosen the quantity most easily accessible experimentally - the 
El amplitudes. It also gives us an opportunity to look for experimental signatures of chaos 
in the Ce atom using the work by Bisson et al ||15||, where over 200 line strengths were 



measured for transitions between a large number of levels within 3.5 eV of the ground state. 
It should be mentioned that there are many other possible atomic systems to search for 
quantum chaos, e.g., in doubly excited states and inner-shell excitation spectra of alkaline- 
earth atoms JIE|-|I8|], or even multiply excited states of light atoms |T9 . 



A. Chaotic many-body states 

Let us now recall briefly what chaotic many-electron atomic eigenstates are. Suppose 
one uses a basis of some single-electron orbitals (e.g., the Hartree-Fock ones) to construct 
many-electron basis states !$&:)• The states can be taken as single-determinant states 
corresponding to certain configurations of a few valence electrons, or constructed from them 
through some coupling scheme to be of definite total angular momentum J. The true atomic 
eigenstates 

l**> = E^°l**>. (E^ )2 = i), a) 



and eigenvalues E^ are obtained by diagonalizing the Hamiltonian matrix Hjk = {<$>j\H\<&k). 
The coefficients C% describe mixing of the basis states by the residual Coulomb interaction. 
In the multi-electron excitation range the number of basis states \<f>k) formed by distribut- 
ing several electrons among a few open orbitals is large. Many of these states are nearly 
degenerate and the mean spacing between the basis state energies Ek = H^k is likely to be 
smaller than the typical value of the off-diagonal matrix element Hjk. In this situation the 



basis states are strongly mixed together p0 | 



Apart from a few lowest levels, each of the eigenstates is a superposition of a large number 
of basis states. Of course, by a simple perturbation theory argument, the mixing must be 
weak for distant basis states (large \Ej — Ek\). The strong mixing takes place within a certain 
energy range \Ej—Ef.\ ^5 T = 2ttV 2 /D, where D is the mean level spacing, V 2 = H 2 k , and T is 
called the spreading width, since it characterizes the spread of the eigenstates to which a given 
basis state contributes noticeably. One can estimate the number of principal components, i.e. 
those that contribute significantly to a given eigenstate (P, as N ~ T/D. The coefficients 
Cjp corresponding to the principal components have typical values \Cj. | ~ l/y/N. Their 
statistics is close to that of independent random variables, and tends towards Gaussian 
when the mixing is strong. In this case even the single-electron orbital occupancies are far 
from integer and only the total angular momentum, parity and the energy itself remain 



good quantum numbers ||11|| . Thus, we can talk about quantum chaos in the system. This 
situation is similar to that in compound nuclei and the corresponding chaotic eigenstates 
can be called atomic compound states. The model configuration interaction calculations 
performed for Ce produced a value of T ~ 2 eV, and demonstrated the existence of a dense 
spectrum of chaotic compound excited states with N > 100 (D ~ 0.01 eV) just few eV from 
the ground state |TT . 
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II. MATRIX ELEMENTS BETWEEN CHAOTIC STATES 



Consider two chaotic many-body states (compound states, for short) that are superpo- 
sitions of large numbers of basis states, = Y^k^i 1^*) an d 1^2} = HjCf^Qj). If the 
expansion coefficients are random, the matrix element of some operator M 



(* 2 \M\* x ) = Y,Cf {^\M\^ k )Cl 

jk 



(1) 



(2) 



is a sum of a large number of almost uncorrelated random items Therefore, one should 
expect that such matrix elements display Gaussian statistics with zero mean. Hence, the 
probability distribution of the matrix elements between compound states can be character- 
ized by their mean squared value alone. 

If M is a single-particle operator, e.g., the electric dipole moment D = Y, a p( a \d\f3)a> a a l3 
(a and (3 are single-particle states), it is convenient to express its matrix elements in terms 
of the matrix elements of the density matrix operator p a p 



(3) 



a/3 



a/3 



where p% ] 



(^ 2 |p Q/3 |^i). 



In 



14J and |l l|] a statistical approach to calculation of mean squared matrix elements 
between compound states has been developed. It is first based on the assumption that 
contributions from different single-particle transitions (3 — > a in the matrix element (|3]) are 
uncorrelated. The mean squared value is then given by 



\(^2\D\^ 1 ) 



\d. 



a3 



2 \o (21) \ 

\Paf3 I 



(4) 



aj3 



where averaging is done over a number of compound states around and/or ^2- The mean 

squared value of the density matrix operator | 2 is expressed in terms of the parameters 
of the compound states 1 and 2 (i.e. their energies and spreading widths), and the average 
occupation numbers of the single-particle states a and (3. 

In a spherically symmetric system where the states 1 and 2 are characterized by their 
total angular momenta Ji^ and projections M\^, the Wigner-Eckhart theorem applies, and 
it is convenient to deal with the reduced matrix elements (\E , 2||-D|| 1 I / i) independent of the 
projections M\^- For example, the mean squared value of the zero-rank reduced density 
matrix operator (Ji = J 2 = J then) is obtained in the following two forms p"Tf : 



J21)0 
rnlj,n'l'j 



d x 5(r 1; r 2 ,A) 
d 2 5(r!,r 2 ,A) 



2J+1 



2.H 
2JH 



(n n ij(l 



2j+l 



'•■n'l'j 
2j + l 



)h 



(5) 



where D 12 are the mean level spacings near the states 1 and 2, n n ij and n n >i>j are the orbital 
occupation numbers, and 5 is a "finite-width 5 function". It depends on the spreading 
widths 1^2 of the compound states and on the energy difference A = uj n 'i'j,nij ~~ 
between the transition frequency for the compound many-electron states E^ 1 ' — and 
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the frequency u n 'i'j,rdj °f the single-particle transition between the orbitals nlj and n'l'j. 
The function 5 has a maximum at A = and describes the energy conservation for the 
compound states. Its width is determined by the spreading widths r 12 . Note that (. . .)i >2 
in Eq. @ denote averaging of the occupation-number factors over the compound states 1 
or 2. Note also that the exact form of the function S(Ti,T2, A) depends of the spreading 
of the compound states over the basis components, i.e. on the "shapes" of the eigenstates. 
In the simplest approximation this spreading is described by the Breit-Wigner formula (see 
numerical studies in [[11]]) and 5 is also a Breit-Wigner profile 

To calculate the mean squared value of the El amplitude we now need a formula for the 



reduced density matrix operator of the first rank. Starting from the definition ] 1 



(21)1 _(_ 1) J 2 -M 2 ( J 2 1 Jl\ V- r _ lV -mf 3 1 f\ (21) (7) 

Pm^l'f-i i; {-M 2 qMj ^} ij {-mqm>) Pnljm > n ' l ' j ' m ' [V 

for q = (linear polarization along the quantization axis) and assuming that transitions 
between different magnetic sublevels m are uncorrelated we can derive a formula for the 
mean square of (0), and then use it to obtain the mean-squared El amplitude, 

o 7 _i_ i 

\(*i\\b\\* 1 )\* = ±±^D % £ |<ny||d||n'J'j'>| a 

nlj,n'l'j' 




2j> + 1 




1-7^7 , (8) 



analogous to the lower formula in Eq. ([|), or an alternative form with Di and 
(27+1 (} ~ 2/+1 )) 011 ^ ne r ight-hand side. The factor 1/3 on the right hand side of Eq. 
(§) is due to the fact that there are three final-state momenta J% — Jx, J\ ± 1 accessible 
from a given J\ by means of a dipole transition. In deriving this expression an additional 
assumption has been made that the occupancies of the nljm and n'l'j'm' states are statis- 
tically independent, and the states with different m within the same nlj shell are equally 
populated. This supposition influences only the "emptiness" factors (l — ^j^j, which are 
close to unity anyway when the number of single-electron states available is much greater 
than the number of active electrons. 

The square of the reduced dipole matrix element S(2, 1) = K^H-DH^i)! 2 is called the 
strength of the line 1 — > 2, so Eq. (||) allows one to estimate mean line strengths for transitions 
involving compound states. 

It is interesting to note that the statistical theory expression (||) satisfies the dipole sum 
rule [^] (in atomic units), 

9 R(2) - 

tt E 9 y _i_ 1 !(* 2 |PII*i)i 2 ^, (9) 

where n is the number of active valence electrons included in the configuration space of the 
problem. To obtain this result one should replace summation over the final states 2 with 
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integration over dE^/D 2 , take into account that J(E^ — E^)5(Ti, T 2 , A)dE^ = u>nij,n'i'j' 
[see Eq. @], neglect the "emptiness" factor (l — |^y) ~ 1 and use Y^ n 'i'j'( n n'i'j')i = n, and 
rely on the single-particle sum rules for the orbitals n'Vf occupied in the initial state ^i, 



III. NUMERICAL RESULTS FOR THE CE ATOM 
A. Energy levels 

Cerium, Z = 58, is the second of the lanthanide atoms. Its electronic structure consists of 
the Xe-like Is 2 . . . 5p 6 core and four valence electrons. The atomic ground state is described 
by the 4/6s 2 5g? configuration with J 71 = 4r [[10 . 



The origin of the extremely complex and dense excitation spectra of the rare-earth atoms 
is the existence of several open orbitals near the ground state, namely 4/, 6s, 5d, and 6p, 
or, in relativistic notation, 4/ 5 / 2 , 4/ 7 / 2 , 6si/ 2 , 5d 3 / 2 , 5d 5 / 2 , 6pi/ 2 , and 6p 3 / 2 . These make 
a total of N s = 32 single-electron states. For Ce with n = 4 valence electrons there are 
about (N s ) n /n\ « 4 x 10 4 possible many-electron states constructed of them. If we allow for 
the two possible parities, about ten possible total angular momenta J, and 2 J + 1 different 
projections (another factor of ten), there will be still hundreds of energy levels within a given 
J 71 manifold. 

In the present work we perform relativistic configuration interaction calculations in the 



Hartree-Fock-Dirac basis analogous to those in |L_TJ. In that work we limited ourselves to 
just 7 nonrelativistic configurations constructed of the 4/, 6s, 5c?, and Qp orbitals, for both 
odd and even states, which produced 260 and 276 states with J n = 4~ and 4 + , respectively. 
To make the results more realistic we have extended the configuration basis set by 9 odd and 
23 even nonrelativistic configurations. Basically, the additional configurations were obtained 
by exciting one of the four electrons of an "old" configuration into the next orbital, e.g., the 
lowest even 4/ 2 6s 2 configuration would produce 4/6s 2 5/, 4/6s 2 7p, 4/ 2 6s7s, and 4/ 2 6s6rf 
configurations. To keep the size of the configuration space reasonable we included only 
the configurations with mean energies within about 10 eV from the Ce ground state. This 
increased the total number of 4~ and 4 + states to 862 and 1433, respectively. Note that 
J = 4 states have been chosen because these manifolds are among the most abundant. 

As a result, the level density p(E) = Y,i ${E — E^) has increased greatly above 4 eV and 
become closer to that observed experimentally. Of course, to be meaningful the level density 
must be averaged over some small energy interval to obtain a smooth function rather than 
a set of spikes. An alternative procedure is to look at the cumulative number of levels 

N(E) = [ E P {E')dE' , (11) 



which we present in Fig. [I] for J n = 4 + states. Each N(E) plot is a staircase of steps of the 
unit height occurring at successive excited state energies. The level density can be easily 
estimated from the slope of the N(E) plot. The experimental data for the 132 even levels 
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with J = 4 known from [|l(]] is shown by the solid-line staircase, and the energies are given 
with respect to either experimental, or calculated ground state energy. They can be com- 
pared with the dashed line that shows N(E) for our earlier small-basis calculation [|11[] (276 



states), and the dotted line for the present calculation (1433 states). The improvement is 
obvious, however, the agreement is not perfect. We believe that the remaining disagreement 
is not due to some missing configurations in the CI calculation, but rather due to an overall 
"softening" of the spectra due to screening of the Coulomb repulsion between the valence 



electrons by the electrons of the core |p2]| . In the CI language this effect is produced by 
the high-energy excitations of the valence electrons into the continuum together with the 
electron excitations from the core. 

Two typical features can be observed in the spectra of complex atoms |J . The first clearly 
seen in Fig. [I], is the rapid increase of the level density p(E) with energy p3]| . Its origin 
is purely combinatorial - the larger the excitation energy, the greater the number of ways 
it can be distributed among a few single-particle excitations. In the independent-particle 
model this dependence is described by the following exponent 0]: 



Pa (E)= Po exp(aJE-E g ), (12) 



where p and a are some constants, and E g is the ground state energy of the sys- 
tem. This dependence also follows from the thermodynamic definition of the temperature, 
T _1 = d{ln[p(E)]} / dE , combined with the estimates of the average number of excited Fermi 
particles, n ex oc T, and that of the excitation energy, E — E g ~ n cx T. The experimental 
spectra of rare-earth atoms and their ions examined in || are in agreement with Eq. (|i~2|). 
Figure [I] shows that the calculated cumulative level number plot is fitted well by 

N(E)= f E p a (E')dE' (13) 



E a 



with p = 0.65 eV -1 , a = 2.55 eV 1//2 , and the "ground state" energy of the 4 + sequence E g 
shifted by 0.25 eV up from the true J n = 4~ ground state of Ce. Thus, Eq. (|T2] ) gives a 
good overall fit of the calculated level density below 6 eV. 

The second feature typical for the spectra of complex many-body systems is level repul- 
sion. It is a basic quantum mechanics fact that two levels with identical quantum numbers 
cannot be degenerate if they are coupled by a non-zero matrix elements - they "repel" 
each other. In quantum chaotic systems this repulsion is characterized by the Wigner level 
spacing distribution 

P(s) = ^e- 2 / 4 , (14) 

where s is the nearest-neighbour level spacing normalized so that s = J sP(s)ds = 1. 
Equation ([14]) shows that the probability of finding small level spacings is indeed vanishingly 
small. As we pointed out in the Introduction, spectral statistics do not tell much about the 
eigenstates of the system. However, Eq. (|i~4]) is still a good test for some possible hidden 
quantum numbers, e.g., the total spin or orbital momentum, which might characterize atomic 
eigenstates besides J n . If these do exist, small level spacings ("degeneracies") will be more 
abundant than predicted by Eq. ([14]). These statistics were checked for many experimental 
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Pj6| JTl|JT6|JT^| and calculated fTTUTTI] complex atomic spectra, as well as for molecular vibronic 



spectra 24 



As seen from Fig. [I] the level density changes significantly for the first 500 levels of the 
calculated spectrum. To analyze the distribution of the corresponding level spacings we use 
the analytical density fit p a {E) to normalize the spacings: 

s n = (E n+ i - E n )p a (E n ). (15) 

Their distribution shown on the inset in Fig. [TJ is in reasonable agreement with the Wigner 
formula. The deviations are probably due to the long-range fluctuations of the level density, 
not accounted for by the simple exponential (|T~2|) . In the previous calculation ||11|| , where 
only the lowest orbitals of each symmetry were included, we also observed the Wigner 
distribution. When orbitals with higher principal quantum numbers become involved (as 
seen from Fig. [I] above 3.5 eV) the spatial extent of the eigenstates increases. This should 
cause a decrease of the residual Coulomb interaction between the electrons. On the other 
hand, the level spacings also become smaller. As a result, the state mixing at these excitation 
energies remains strong, which is confirmed by the agreement with the Wigner distribution, 
and the eigenstates are chaotic. Our estimate of the number of principal components N 
shows that it becomes even greater as the energy increases, in accord with the estimate 
N ~ Y/D ~ 300 (r ~ 1 eV, and the mean level spacing D « 0.003 eV at E w 6 eV). 



B. Dipole matrix elements 

In Sec. H we explained that matrix elements involving chaotic compound states should 
have Gaussian statistics, and the mean squared value of the matrix elements could be esti- 
mated in terms of some average characteristics of the compound states. In this section we 
concentrate on the dipole matrix elements (.El amplitudes) d^ = (tyf \\D\\tyf, ) between 
the 14 lowest states with J n = 4~ and 80 consecutive 4 + states obtained numerically in our 
CI calculations of Ce. We have chosen this energy region to cover the range explored in 



the experiment [15], where absolute values were derived for 228 of the most intense lines of 
neutral Ce between 10706 and 22184 cm -1 . 

Of course, low-lying atomic states, e.g., the ground state, have well-defined configuration 
composition and are not chaotic, hence, the El amplitudes between them should not be 
distributed in any particular statistical way. However, the matrix elements (0) will become 
random (and close to Gaussian) as soon as at least one of the states involved, the initial 
or the final, moves into the compound-state energy range and becomes a superposition of 
many random components. Besides that, the mean squared value of the matrix element 
is expected to show some smooth secular variation with the energy of the states involved. 
For these reasons we skip the first 20 states with J n = 4 + and analyze the statistics of the 
14 x 80 = 1120 El amplitudes for the following 80 even states by grouping them in bunches of 
twenty: 21-40, 41-60, 61-80, and 81-100, which correspond to the mean excitation energies 
of 2.49, 2.95, 3.40, and 3.70 eV above the atomic ground state (the mean energy of the 
lowest 14 odd states is 0.68 eV). Thus, each plate in Fig. |2] shows the distribution of the 
280 reduced dipole matrix elements together with their root-mean-square (r.m.s.) value. 



Also shown in Fig. ^|are the Gaussian distributions g(d) = exp(—d 2 /2dl)/J2iTdl, where the 
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root-mean-square parameter do has been adjusted to minimize x 2 around the center of the 
histogram. The values of do and x 2 are given in Table |. 

Two effects can be seen in Fig. |2|. First, the distributions of the matrix elements are 
indeed close to Gaussian. Second, the width of the distributions (the mean squared value of 
the matrix elements) varies with the energy of the even states. It is mostly this effect that 
is responsible for the visible discrepancies between the histograms and the Gaussian fits. To 
eliminate it we can use a running average procedure to normalize the amplitudes: 

where {d 2 )}^ 2 is the r.m.s. value over the 14 odd states, calculated for every even state i. 
Figure § confirms that the 1120 normalized El amplitudes for the 21-100 even states are 
distributed according to the normal law. The inset shows the dependence of the r.m.s. El 
amplitude (d 2 )]^ 2 on the energy of the even state E^\ Fluctuations aside, it is in agreement 
with the r.m.s. values calculated from the statistical theory, Eq. (^), at the energies of the 
30th, 50th, 70th and 90th even states. The numerical values of the r.m.s. El amplitudes 
are listed in Table |. 

Note that we have chosen Eq. (pp with 1 standing for the odd states and 2 for the even 
ones. In our numerical example we consider the dependence of the r.m.s. El amplitude on 
the energy of the even states, and keep the odd states the same. Therefore, as in Eq. 
we only need to know the average occupation numbers for the lowest 14 odd states, and 
the result depends on the final even state via its energy E®, mean level spacing D2 and 



spreading width T 2 . As we saw in our previous calculations QTTJ , the even states of Ce with 



J = 4 become very much chaotic at excitation energies of just 2 eV, i.e. from the 20th level 



up. Also, as earlier in we use average configuration energies rather than single-particle 
Hartree-Fock energies to determine the transition frequencies to n >i>j, n ij needed for calculation 
of A in Eq. (^). The ground state of Ce is described as 4f6s 2 5d, however, the dominant 
configuration among the 14 lowest odd states is 4f6s5d 2 , and we used it to calculate the 
transition energies. For example, the energy of the 6s — 6p transition ujq p ^ s was determined 
as the difference between the average energies of the 4f5d 2 6p and 4fQs5d 2 configurations. 
Physically, this corresponds to choosing a particular mean field close to that of the low-lying 
odd states of Ce for calculation of the transition energies. It should be mentioned though 
that the results obtained with the Hartree-Fock frequencies uj n ij, n 'i'j' — e nij ~ e n'i'j' were not 
too different. 

Gaussian statistics of the dipole matrix elements result in the Porter-Thomas (PT) dis- 
tribution of the line strengths S(i, k) = d 2 k 

f(S) = ^^e W (-£=), (17) 



V2^ V 2S 

where S is the mean line strength. Divergence of this function at small S means that if the El 
amplitudes are Gaussian, there should be many weak lines in the spectrum. Earlier evidence 
of the PT statistics of line strengths can be found in calculations of dipole excitations in 



complex atoms ||25|| , and transitions between the vibronic levels in molecules measured in 
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IV. ANALYSIS OF EXPERIMENTAL DATA 



In [13] absolute values of gA were obtained for 228 of the most intense observed lines 
between 10706 and 22184 cm -1 in Ce. It is interesting to analyze these data to see whether 
they support our theoretical and numerical considerations. 



The values of gA listed in |0J are defined as gA = (2Jk + l)Aki, where 



Aki = ff^ \(i\\D\\k)\ 2 (18) 

3ftc 3 (2J fc + l) h 11 11 71 K J 



is the El transition rate from the upper level k into the lower level i ||26|| . We use the 
experimental values of gA, Jk and transition frequencies lo^ to extract values of the line 
strengths 

s(i,k) = \m\\k)\ 2 = gA-f¥j.. (19) 

In Fig. |] the probability distribution of the 228 experimental line strengths is shown. 
Compared to the expected PT formula ([H]), there is a clear lack of small line strengths. 
Nevertheless, the decreasing part of the histogram can be fitted well by a PT distribution 
with an additional normalization factor A, 

= A^i-SftS) , (20) 
V2irSS 

shown in Fig. [| by a solid line for A = 2.07 and S = 3.3 a.u. that minimize \ 2 f° r the 22 
bins with S > 3 a.u. 

It would be tempting to say that the excellent agreement between the PT curve and the 
histogram is a confirmation of the Gaussian statistics of the El amplitudes in Ce. The value 
of A would then indicate that about one half of all lines are missing in the experimental 
data. However, the value of S — 3.3 a.u. corresponds to the r.m.s. El amplitude of 1.8, 
which is more than 2 times greater than our numerical results in Fig. ||] and Fig. || (inset) 
and in Table |. On the other hand, the experimentally observed 228 lines include transitions 
between levels with various total angular momenta between J = 1 and 8 (| Jj — J^\ < 1, of 
course) whereas we have about 500 hundred lines with just Ji = Jk = 4 in our calculation 
in the analogous energy range. This means that in |T5J only the strongest 10% or less of all 
lines have in fact been measured. The very suggestive agreement with the PT distribution 
in Fig. f| should then be considered as merely fortuitous. 

It is worth noting that in experiment the lines are selected by their intensities propor- 
tional to gA, rather then strengths. Hence, even lines with large strengths can be omitted 
if their frequencies are small. Let us look at the simplest model of this effect and see how 
it influences the observed strengths distribution. Assume that transitions in a certain fre- 
quency range < uj < u max are studied, and different values within this interval are equally 
probable. The observed intensities of the lines are proportional to u 3 S. If we assume that 
there is a minimal threshold intensity that can be registered, the original PT distribution of 
strengths would be modified as follows: 
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, S < So 

A exp ( = ) 



1 



t) 



1/3' 



S ^> So 



(21) 



w mitt , and A is the normal- 



where So is the minimal strength that can be observed at uo 
ization factor. As seen from Fig. |], Eq. fl2~l|) also gives a very good fit of the experimental 
data with S = 2.4, So = 1.12 and A = 6.22 [2~7|. Note, however, that the new value of the 
mean line strength is 1.5 times smaller than the one we had from the pure PT fit. There- 
fore, the assumptions used in our processing of the experimental data affect the estimates 
of the experimental r.m.s. El amplitudes, and we should not be too concerned about the 
apparent disagreement with our numerical calculations. Besides that, extraction of absolute 
line strengths from the experimental data is not free from uncertainties estimated in [15| at 
10-20%. 

For 30 transitions in Ce the gA values were obtained more accurately from branching 
ratios and delayed photoionization measurements of lifetimes ( ||15|| , Table 2). When we 
look at the statistics of the corresponding line strengths (Fig. [|, inset) and compare it with 
the PT distribution (|T7|), a value of S = 2.15 is obtained, much smaller than the estimates 
of S from the statistics of the 228 lines. Thus, it appears that to make firm conclusions 
about Gaussian statistics of the El amplitudes a much more thorough experimental survey 
is needed. On the other hand, even relative measurements of a large number of line strengths 



could be very valuable for examining these statistics p4|. 



V. CONCLUSIONS 



In this work we have extended the configuration interaction approach of |TT| to calculate 
large numbers of eigenstates in Ce. In agreement with our earlier studies the energy level 
statistics indicate that the simple configurational basis states are strongly mixed together 
by the residual electron interaction, and the only good quantum numbers in the spectrum 
are parity and the total angular momentum. The total orbital momentum L and spin S are 
not conserved due to the spin-orbit interaction, whose effect is dynamically enhanced, just 
as that of any other perturbation in a chaotic many-body system ||1 1|| . 

The strong configuration mixing makes multielectron atomic eigenstates chaotic. This 
in turn results in a Gaussian statistics of the matrix elements for chaotic atomic eigenstates 
(compound states). This understanding is fully confirmed by our numerical calculations of 
the 1120 El amplitudes between the 14 lowest J n = 4~ states and 80 J n = 4 + states above 
2 eV. It is important that the parameter of the Gaussian, the r.m.s. El amplitude, varies 
slowly with the excitation energy. This effect should be taken into account when analyzing 
the statistics of the matrix elements. 

We also show that a statistical theory can be used to estimate mean squared matrix 
elements involving compound states. It enables one to express the answer in terms of the 
single-particle matrix elements and occupation numbers, and parameters of the compound 
states, namely the number of principal components and the spreading width. This approach 
has already been applied to calculation of matrix elements between compound states in 
nuclei jnj. It could be useful in various other many-body systems, e.g., atomic clusters 



or quantum dots, where direct diagonalization of the Hamiltonian matrix is not feasible 
because of a huge size of the Hilbert space of the problem. 
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An attempt has been made to analyze existing experimental data for the line strengths 
in Ce [15]. It appears that the statistics of the measured line strengths is compatible with 
the Porter-Thomas distribution, with allowance for the missing weak lines. However, the 
discrepancy between the calculated r.m.s. El amplitudes and those inferred from the ex- 
perimental data does not allow us to say that the existence of quantum chaos in the Ce 
eigenstates has been confirmed experimentally. To make this statement one would have to 
do a much more complete survey and statistical analysis of the line strengths in the Ce 
spectrum. 

On the other hand, this means that a comparison between the experimental and the- 
oretical line strengths in Ce is not yet possible, even at the level of their mean values. 
Theoretically, to calculate precisely the dipole matrix elements between particular levels in 
the compound-state energy range of complex atoms like Ce looks a prohibitively difficult 
problem. Experimentally, identification of specific lines in enormously complicated spectra 
is also a very difficult task. However, we would like to suggest that extraction of mean 
characteristics from the experiment and comparison with the corresponding theoretical esti- 
mates is a meaningful way of exploring such complex systems. As a result, one might hope 
to get a deeper insight into the existence of quantum chaos in many-body systems on the 
whole, and in complex open-shell atoms, in particular. 
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TABLES 



TABLE I. Root-mean-square El amplitudes for transitions between the 14 J 71 " = 4 and 80 
J 71 " = 4+ states in Ce. 



even levels 


r.m.s. 


El amplitudes (a.u.) 
from Eq. (||) b 


d c 


n d 


X 2 (n-1) 


21-40 


0.853 


0.813 


0.729 


21 


24.4 


41-60 


0.891 


0.746 


0.824 


21 


23.9 


61-80 


0.736 


0.674 


0.671 


17 


36.6 


81-100 


0.627 


0.566 


0.632 


17 


24.7 



a Obtained directly from the CI calculation. 

b Calculated from the statistical theory, Sec. [fi]. 

c Values that minimize x 2 f° r the Gaussian fits shown in Fig. |2[ 

d Number of bins around the centre of the histogram used for calculation of x 2 ■ 
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FIGURES 



FIG. 1. Energy spectra and level statistics of the J 77 = 4 + states in Ce. Dashed line shows the 
cumulative number of states N(E) for the calculation with 276 basis states dotted line is the 
present calculation with 1433 basis states; thick solid line is N(E) for 132 experimental levels from 
[10]. Thin solid line is the cumulative level number corresponding to the independent-particle fit 
(|13|). Shown on the inset are the statistics of the normalized level spacings s for the lowest 500 
levels, compared with the Wigner distribution (14). 



FIG. 2. Probability distributions of the El amplitudes in Ce for transitions between the 14 
lowest 4" states and groups of 20 states with J n = 4+: (a) 21-40, (b) 41-60, (c) 61-80, and (d) 
81-100. R.m.s. values of the amplitudes are shown next to the histograms. Smooth curves are 
Gaussian fits that minimize x 2 for 21 to 17 central bins of the histograms (see table |). 



FIG. 3. Probability distributions of the normalized El amplitudes in Ce for transitions between 
the 14 lowest 4 _ states and 21-100 states with J 71 " = 4 + , compared to the normal distribution (solid 
line). The inset shows the dependence of the r.m.s. El amplitude averaged over the 14 odd states 
on the energy of the even state (thin solid line). Solid circles connected by a thick solid line are 
the r.m.s. values of the El amplitude obtained from the statistical theory, Eq. (g), at the energies 
of the 30th, 50th, 70th, and 90th even states, and open circles are values from the CI calculation 
(see table |) . 



FIG. 4. Comparison of the line strengths measured in Ce by Bisson et al [15| with the 
Porter-Thomas and modified Porter-Thomas distributions. Solid line is the PT distribution (20) 
with A = 2.07 and S = 3.3 a.u., and dashed line is the modified PT distribution ( |2l| ) with A = 6.22, 
Sq = 1.12 and S = 2.4 a.u. Shown on the inset is the probability distribution of the 30 lines mea- 
sured from branching ratios and delayed photoionization, fitted by a PT distribution with S = 2.15 
a.u. 
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